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ARE THERE FIELD-FREE GAPS NEAR r= 1 IN SUNSPOT PENUMBRAE ? 
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ABSTRACT 

The vertical stratification of the magnetic field strength in sunspot penumbrae is investigated by means 
of spectropolarimetric observations at high spatial resolution from the Hinode spacecraft. Assuming 
that the magnetic field changes linearly with optical depth we find that, in those regions where the 
magnetic field is more inclined and the Evershed flow is strongest (penumbral intraspines), the magnetic 
field can either increase or decrease with depth. Allowing more degrees of freedom to the magnetic 
field stratification reveals that the magnetic field initially decreases from logTs = —3 until logrs ~ 
— 1.0, but increases again below that. The presence of strong magnetic fields near the continuum is at 
odds with the existence of regions void of magnetic fields at, or right below, the T5 = 1 level in the 
penumbra. However, they are compatible with the presence of a horizontal flux-tube-like fleld embedded 
in a magnetic atmosphere. 

Subject headings: Sun: sunspots - Sun: magnetic fields - Sun: polarimetry 



1. INTRODUCTION 

It is now widely accepted that the horizontal structure of 
the sunspot penumbra is composed of two magnetic com- 
ponents (Solanki 2003; Bellot Rubio 2003). One of them 
possesses a somewhat inclined (~ 40 — 50° with respect 
to the vertical direction to the solar surface) and strong 
2000 G) magnetic field, whereas the other is charac- 
terized by a weaker and more horizontal one (Lites et al. 
1993; Riiedi et al. 1998; Bellot Rubio et al. 2004; Bor- 
rero et al. 2004, 2005). Traditionally, these two mag- 
netic components have been identified with a horizontal 
flux tube, that carries the Evershed flow, and is embed- 
ded in a more vertical background magnetic field: un- 
combed model (Solanki & Montavon 1993; Schlichenmaier 
et al. 1998; Borrero 2007). Recently, this view has been 
challenged by Spruit & Scharmer (2006) and Scharmer 
& Spruit (2006), who propose instead that the penum- 
bra is formed by magnetic field-free plumes (connected to 
the underlying convection zone) that pierce the penumbral 
magnetic field from beneath. This is the so-called gappy 
penumbral model. 

So long as these two different magnetic structures 
(weak/horizontal and strong/vertical) have remained spa- 
tially (horizontally) unresolved, distinguishing between 
the uncombed and gappy penumbral scenarios has not 
been possible. However, with the new spectropolarime- 
ter on board of the Japanese spacecraft Hinode (Kosugi 
et al. 2007; Shimizu et al. 2007) it is now possible to ob- 
tain high spatial resolution (~ 0.32") observations of the 
sunspot penumbra. This could be sufficient to distinguish 
between the uncombed and gappy models, since they both 
postulate the existence of fiux tubes or field-free gaps that 
are about 200-300 km in diameter (Martinez Pillet 2000; 
Spruit & Scharmer 2006). This feature is particularly in- 
teresting, because the uncombed and gappy models predict 



very different vertical stratifications in the magnetic field 
strength across the weak/horizontal magnetic field com- 
ponent: which is identified with an embedded flux tube 
in the uncombed model, but with a field-free gap in the 
gappy model. In the latter, the magnetic field decreases 
monotonically with depth, whereas the former possesses 
a magnetic field that decreases with depth only initially, 
since once the boundary of the fiux tube is reached, the 
magnetic field can either decrease or increase depending 
upon the strength of the magnetic field inside the tube. 

In this paper we will focus on obtaining the vertical 
stratification of the magnetic field for penumbral filaments 
(where the magnetic field is more horizontal and weaker) 
using high spatial resolution spectropolarimetric observa- 
tions from Hinode, in order to establish which penumbral 
model is more realistic. The observations are described 
in Section 2. Section 3 describes our data analysis and 
results from our inversion technique. Section 4 compares 
our findings with the predictions made by the uncombed 
and gappy penumbral models. In Section 5 we make a 
thourough investigation of the effects of the scattered light. 
Finally, Section 6 summarizes our findings. 

2. OBSERVATIONS 

On May 3rd 2007, between 10:15 and 11:40 am UT, 
the active region AR 10953 was mapped using the spec- 
tropolarimeter of the Solar Optical Telescope on-board of 
the Hinode spacecraft (Lites et al. 2001). The active re- 
gion was located at a heliocentric angle oi 9 ~ 19.2°. It 
was scanned in a thousand steps, with a step width of 
0.148" and a slit width of 0.158". The spectropolarimeter 
recorded the full Stokes vector (/, Q, U and V) of the pair 
of neutral iron lines at 630 nm with a spectral sampling 
of 21.53 mA. The integration time was 4.8 seconds, result- 
ing in an approximate noise level of 1.2 x 10"'^ (in units 
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of the normalized continuum intensity). In the absence of 
the telluric oxygen lines we proceeded with two different 
wavelength calibration methods that were cross-checked 
for consistency. The first method was obtained by match- 
ing the average quiet Sun profile with the FTS spectrum, 
whereas the second calibration assumes that the average 
umbral profile exhibits no velocities. 

A map of the continuum intensity at 630 nm of the 
scanned region is shown in Figure 1. The white arrow 
indicates the direction of the center of the solar disk. The 
penumbra on the center side is heavily distorted and there- 
fore left out from our analysis. On the limb side the 
penumbra is more uniform, with radially aligned filaments. 
The region enclosed by the white rectangle has been chosen 
for our study. This sunspot has negative polarity (mag- 
netic field in the umbra points towards the solar interior), 
however the results presented hereafter are shown, in or- 
der to facilitate the interpretation, as if the sunspot had 
positive polarity. 
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Fig. 1. — Continuum intensity map at 630 nm of AR 10953. This 
sunspot was observed using Hinode's spectropolarimeter on the 3rd 
of May, 2007 at an heliocentric angle oi 6 = 19° . The white arrow 
points towards the center of the solar disk. The white rectangle lim- 
its the region chosen for our study. It lies on the limb-side around 
the line-of-symmetry. 

3. DATA ANALYSIS AND RESULTS 

We have applied the SIR inversion code (Ruiz Cobo & 
Del Toro Iniesta 1992) to our spectropolarimetric obser- 
vations to retrieve the physical properties of the solar at- 
mosphere. This code allows all relevant physical parame- 
ters to be a generic function of the optical depth: B{t), 
7(t), 4>{t), T4os('''), etc. In addition, a depth-dependent 
temperature stratification T{t) models the atmospheric 
thermodynamics under local thermodynamic equilibrium 
(LTE) conditions. SIR retrieves the values of the param- 
eters at a number of optical depth points called nodes. 
The final stratification is obtained by interpolating splines 
across those nodes. Note however, that SIR employs equiv- 
alent response functions (Del Toro Iniesta 2003), which en- 
sures sensitivity to the atmospheric layers located between 
nodes. Each node represents a free parameter in the in- 
version. In our investigation we will employ increasingly 
complex models (i.e.: more free parameters) according to 
the amount of information we hope to extract from the 
profiles. 

Given the high spatial resolution of Hinode's observa- 
tions, we will consider only one magnetic component. A 



non-magnetic component is also considered to account for 
the scattered light. In this section, the scattered light 
profile is obtained by averaging the intensity profiles of 
those pixels with polarization signals below the noise level 
(quiet Sun granulation around the sunspot). The same 
scattered light profile is used in the inversion of all pixels. 
In Sect. 5 we make a thourough analysis of the effects that 
different treatments for the scattered light have on our 
results. Note that using one single magnetic component 
is equivalent to assuming that the penumbral structure is 
horizontally resolved. This is clearly not the case if we 
look into continuum images at even higher spatial reso- 
lution (Scharmer et al. 2002). However, our assumption 
would still be valid if the (unresolved) variations of the 
magnetic field inside the weak/horizontal magnetic com- 
ponent are much smaller than the differences between the 
weak/horizontal and strong/ vertical components. Since 
the former remain unresolved, we cannot assess the valid- 
ity of this assumption. This question should be addressed 
as better spectropolarimetric observations become avail- 
able. 

3.1. 1-node inversion and intraspine selection 

In order to locate the intraspinal pixels we have carried 
out a first inversion where all physical parameters, with 
the exception of the temperature, are constant with optical 
depth. We therefore have one single node for -B(t), 7(t), 
4'{t)^ Mos('''). To account for unresolved velocity fields, 
we also consider depth-independent micro and macrotur- 
bulent velocities: Vmic and Knac- Another free parame- 
ter, ttqs, represents the fraction of the observed intensity. 
Stokes /, that corresponds to scattered light. Finally, three 
nodes are given to the temperature T{t). In total, this first 
inversion has 10 free parameters. Since B, 7, and Vios are 
constant with optical depth, the retrieved values indicate 
some kind of average over the region where the spectral 
lines are formed. Westendorp Plaza et al. (1998,2001) 
studied this issue in detail and found that the largest con- 
tribution for this pair of Fe I lines (Sect. 2) comes from 

l0gT5 ~ -1.5. 

Figure 2 displays the resulting values for the line-of-sight 
velocity and magnetic field vector in the selected box in 
Fig. 1. Regions of weak, B < 1300 G, and highly inclined, 
7 > 80°, magnetic field can be clearly distinguished in 
this figure. They are also characterized by the presence of 
large red-shifted velocities (Evershed flow). These are the 
so-called penumbral intraspines, and therefore the most 
likely locations where field-free gaps or horizontal fltix- 
tubes can be found. Also visible are structures charac- 
terized by a stronger and more vertical magnetic field, as 
well as by a strongly reduced Evershed flow. These are 
usually referred to as spines. Spines and intraspines are 
also seen at moderate (~ 1") spatial resolution (Lites et 
al. 1993; Stanchfield et al. 1997; Mathew et al. 2003) but 
the associated changes in their properties (field strength, 
inclination, etc) are larger if observed at high spatial reso- 
lution (Bello Gonzalez et al. 2005; Langhans et al. 2005). 
Bellot Rubio et al. (2004) interprets this result as a conse- 
quence of these structures not being spatially resolved at 
1" resolution. As demonstrated by Borrero et al. (2008) 
they are indeed horizontally resolved in Hinode observa- 
tions (0.32"). 
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Fig. 2. — Results from the inversion of the region hmited by the square box in Fig. 1. The inversion was performed assuming that S, 
7, <j> and Vios are constant with optical depth. The magnetic field strength is displayed in the lower-left panel, inclination (upper-right), and 
line-of-sight velocity (upper- left). The white and black dots correspond to those pixels where the location of horizontal flux tubes or fleld-free 
gaps are suspected (see text for details). There are 7520 of them: 39 % of all penumbral pixels in this figure. 
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Fig. 3. — Temperature (upper-left), line-of-sight velocity (upper-right), magnetic field strength (lower-left), and magnetic field zenith 
angle (lower-right) as a function of the optical depth logrs, obtained from the 2-node inversion of the intraspinal pixels selected in Fig. 2. 
Approximate height scale, Z, computed assuming vertical hydrostatic equilibrium, is also indicated. 
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In Figure 2 wc also indicate with black and white dots 
a large number (total of 7520) of intraspinal pixels. They 
have been found as those satisfying the following condi- 
tions: B > 700 G, and V|os e [1.0,3.0] km s'^ Since 
the main difference between spines and intraspines is the 
presence of a strong Evershed flow, we use Vios to distin- 
guish among them. However, we do not consider the few 
pixels where Vios > 3 km s~^, since they usually present 
extremely abnormal Stokes V profiles, usually a sign of the 
existence of horizontally unresolved structure. We do not 
constrain the values of the magnetic field inclination and 
strength (here we use only a lower limit to avoid taking 
pixels outside the visible boundary of the sunspot) be- 
cause the magnetic properties of the spines in the outer 
penumbra are very similar to those of the intraspines in 
the inner penumbra. The final selected pixels represent 
about 39 % of all penumbral pixels in Figure 2. Note that 
they are mostly located in the middle and outer penum- 
bra: r/Rs > 0.5 (i?s being the sunspot radius; umbral- 
penumbral boundary is located at r/Rg > 0.25). Note also 
that, even though we have not constrained the values of 
the magnetic field strength an inclination, all intraspinal 
pixels are located in regions where the magnetic field is 
highly inclined and weak. 

3.2. 2-node inversion of individual profiles 

In order to investigate the depth variation of the phys- 
ical parameters in intraspines, we performed a renewed 
inversion of the pixels selected in Fig. 2, where we now 

allow for two nodes in -B(t), 7(r), 0(r), V\os(t) (linear 
variations with optical depth). The total number of free 
parameters is now 14. Results from this new inversion are 
presented in Figure 3: T(r) (upper-left), V|os(t) (upper- 
right), B{t) (lower-left), and 7(t) (lower-right). All in- 
verted pixels display similar stratifications of Vios(t) and 
7(r): both increase monotonically with optical depth: 
dVios/dr, dj/dr > 0. The magnetic field strength B{t), 
however, can either increase (in 66 % of inverted pixels: 
4971) or decrease (34 %; 2249 pixels) with optical depth. 
In either case, the retrieved gradient is relatively small: 
\dB/dz\ < 1.5 Gauss km~^. An important feature to note 
is that pixels displaying a decreasing magnetic field to- 
wards deeper layers, dB/dr < 0, are mostly located in the 
inner penumbra, whereas pixels showing dB/dr > are 
mostly found in the outer penumbra (see Figure 4). 
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Fig. 4. — Percentage of the total selected pixels that, at each 
radial distance from sunspot center, show a magnetic field that in- 
creases {dB/dr > 0; solid line) or decreases {dB/dr > 0; dashed 
line) towards the solar interior. 

3.3. ^-node inversion of individual profiles 

We now perform a more complex inversion of the same 
pixels as in Sect. 3.2. In this case we allow for 4 nodes 
in B{r), j{r), (f>{T), and Vios(T). These nodes are located 
at optical depth positions: logrs = [—3.2, —1.8, —0.4, 1]. 
In total, this new inversion has 22 free parameters. Fig- 
ure 5 shows the results from the 4-node inversion of the 
7250 intraspinal pixels selected in Sect. 3.1. The stratifi- 
cations are very similar to those already obtained through 
the 2-node inversion (see Fig. 3). The larger scatter (pixel- 
to-pixel variations) in the 4-node inversion is due to the 
larger amount of free parameters, which are more weakly 
constrained by the observations. 

Since now we allow for 4 nodes to the stratification of 
the magnetic field strength it is not easy to classify our re- 
sults between those where the magnetic field increases or 
decreases with optical depth. To showcase the differences 
between the possible stratifications we have taken sepa- 
rately those pixels where, in the 2-node inversion, showed 
dB/dr < (family 1) or dB/dr > (family 2) and ob- 
tained the averaged stratification for the 2 and 4-node in- 
version. Results for family 1 and 2 are presented in Figures 
6 and 7 respectively. 

Interestingly, the magnetic field strength, that in the 2- 
node inversion showed different gradients, now shows an 
initial decrease up to logTs ~ — 1 (approximately 100 km 
above the continuum level), where it starts to increase 
again towards dc;ep layers. This happens for both fami- 
lies of magnetic structures, and thus could indicate that 
they are indeed closely related. A closer look reveals that 
both families possess a similar magnetic field strength in 
deep layers: B{r = 1) ~ 1250 G, but slightly different 
higher up: B{r = 10"^) ~ 1200 G (family 1; Fig. 6) and 
B{r = 10-3) ~ 900 G (family 2; Fig. 7). This effect 
explains why the 2-node inversion (Sect. 3.2) retrieves dif- 
ferent overall gradients for the magnetic field strength: it 
is due to a large variation in the magnetic field at around 
Ts ~ 10"^ since the magnetic field deeper down is basically 
the same in both cases. 

To further confirm these results we have repeated our 
4-node inversion with the nodes located at slightly dif- 
ferent positions. SIR always places 2 nodes at the up- 
permost and deepest r-locations of the discretized atmo- 
sphere, while spreading the rest equidistantly in between. 
Therefore, to keep the same number of nodes and, at the 
same time, change their r-positions we must change the 
initial and last r-points of the atmosphere. In our first 
set of inversions the atmosphere is discretized between 
logT5 = [—3.2,1]. Changing this to logT5 = [—3.0,1.2] 
and logTs — [—3.7,0.5] would position the 4 nodes at 
[-3.0,-1.6,-0.2,1.2] and [-3.7,-2.3,-0.9,0.5], respec- 
tively. We have inverted all pixels again in these two cases 
and confirmed that our results (Fig. 5,6 and 7) do not 
change. This is due to the use that SIR makes of equiva- 
lent response functions (see Section 3). 

4. DISCUSSION 
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Fig. 6. — Temperature (top left), line-of-sight velocity (top right), magnetic field strength (lower left) , and magnetic field inclination (lower 
right) as a function of the optical depth. Red indicates the average stratification obtained from the individual 2-node inversion of the 2549 
profiles belonging to family 1: dB/dr < (taken from Fig. 3). Green shows the average stratification obtained from the individual 4- node 
inversion of the Stokes vector of the same pixels. 




If we consider a ray passing through the center of an 
intraspine, the gappy and uncombed penumbral models 
predict a very similar stratification of the magnetic field 
strength above the field-free gap or flux tube, but very 
different ones inside them. Figure 8 illustrates some pos- 
sible stratifications predicted by these two models, where 
the upper boundary of the field-free gap and flux tube is 
located at z = 0. Above z — they both share the same 
stratification for the surrounding magnetic atmosphere. 
Here we present two examples, one where the surround- 
ing magnetic field is weak: i?surr — 1000 (dashed line; 
meant to represent the outer penumbra, r/Rs = 0.8), and 
another case where the surrounding field is stronger (solid 
line): Bsurv = 1500 (meant to represent to inner penum- 
bra, r/Rs = 0.4). Note that the magnetic field strength 
in the surrounding atmosphere decreases towards deeper 
layers. This is due to the fact that the vertical component 
of the surrounding field must vanish (or nearly vanish in 
the case of a cusp-shaped boundary) at the flux tube's or 
gap's boundary. These two examples are actual solutions 
of analytical models (Fig. 5 in Spruit & Scharmer 2006; 
Fig. 3 in Scharmer & Spruit 2006; see also Eqs. 33-34 in 
Borrero 2007). Below the boundary of the flux tube or 
field- free gap, z = 0, both models predict a very different 
situation. In the case of the gappy penumbra this region 
is void of magnetic fields: i?gap — (hollow circles). In 
contrast, the uncombed model assumes the existence of a 
flux tube where the magnetic field is strong Btube = 1250 
G (filled circles). 

If we compare Fig. 8 with our 2-node inversion (Fig. 3) 
of intraspinal pixels we find that, on the one hand, the 
gappy penumbral model can only explain the slowly de- 



creasing magnetic field, observed for 34 % of intraspinal 
pixels (family 1), if the rs = 1 level is formed above the 
gap's boundary, otherwise a much more sudden drop would 
be observed (hollow circles in Fig. 8). On the other hand, 
this model does not offer any explanation for the 66 % of 
the intraspinal pixels that present an increasing magnetic 
field strength towards deeper layers (family 2). However, 
the uncombed penumbral model can explain both observed 
situations. It all depends on the strength of the flux tube's 
magnetic field as compared to the magnetic field high 
above it: -Bsurr versus -Btubc- A magnetic field that de- 
creases smoothly towards the interior of the photosphere 
can be explained by a flux tube (of any field strength) 
whose upper boundary layer lies below T5 = 1. If the up- 
per boundary is above T5 = 1, it can also be explained with 
a magnetic field inside the flux tube that is weaker than the 
magnetic field a few hundred kilometers above (solid lines 
plus filled circles in Fig. 8). In addition, a magnetic field 
that increases towards the interior of the photosphere is 
compatible with a flux tube with an upper boundary layer 
above rs = 1, and with a stronger magnetic field than the 
one above (dashed line and filled circles in Fig. 8). 
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Fk;. 8. — Vertical variation of the magnetic field strength across 
the center of a field-free gap (hollow circles) according to the gappy 
pcnumbral model. Same for a flux tube with a magnetic field 
strength of f250 G (filled circles). Note that both models share 
the same stratification above the tube's or gap's boundary {z > 0). 
The solid line represents a situation where the external field is rather 
strong (inner penumbra), while the dashed line corresponds to the 
outer penumbra (weak external field). Also, note that the rg = 1 
level can be shifted horizontally such that the continuum level can 
be formed above the gap/flux-tube or inside them. 

A more complex (4-nodes) inversion of intraspinal pro- 
files indicates that, what appeared as two different families 
of structures using a 2-node inversion, are likely to cor- 
respond to one single kind of magnetic structure, where 
the magnetic field exhibits an initial decrease between 
logTs € [—3, —1.0], but increases between logr^ G [—1-0, 0] 
(see Figs. 6-7). While the gappy model offers no expla- 
nation for this effect, it can indeed be explained by the 
uncombed pcnumbral model, by means of a magnetic field 
whose strength decreases initially but increases once the 
line-of-sight crosses the flux tube's boundary (see Fig. 8). 
Furthermore, although intraspinal families 1 and 2 ap- 
pear to be the equivalent in the 4-node inversion, they 
still present a subtle yet important difference: family 1 
(more commonly found in the inner penumbra; see Fig. 4) 
displays a much stronger initial decrease as compared to 
family 2, which is usually found in the outer penumbra 
(compare lower- left panels in Figs. 6 and 7). 

This can be explained, in terms of the uncombed model 
if, at small-intermediate radial distances, the horizontal 
flux tube possesses a weaker magnetic field than the field 
in the atmosphere in which it is embedded: -Btubc < -Bsurr 
at r/Rs small (compare solid line plus filled circles in Fig. 8 
with green solid in Fig. 6). As we move towards larger ra- 
dial distances, and assuming that the magnetic field inside 
the flux tube remains constant, the surrounding magnetic 
field weakens and falls below the flux tube's field strength: 
-Btubc > -Bsurr at r / Rg large (compare dashed line plus 
filled circles in Fig. 8 with solid green in Fig. 7). Note 
that the assumption that the magnetic field in the flux 
tube remains constant is in agreement with a surrounding 
magnetic field whose strength decays much more rapidly 
towards the outer penumbra than inside the flux tube (see 
Fig. 4 in Borrero et al. 2004; Fig. 6 in Borrero et al. 2005 
and Fig. 4 in Borrero et al. 2006). Here we find that 
this known feature of the pcnumbral intraspines helps to 
explain, within the frame of the uncombed model, differ- 
ences in the stratification in the magnetic field strength 
across instrapines at different radial distances, as deduced 



from high resolution spectropolarimetric observations. 

5. SCATTERED LIGHT CONSIDERATIONS 

One of the most critical issues in the inversion of spec- 
tropolarimetric data is the treatment of the scattered 
light. In order to properly model its contribution, detailed 
measurements of the telescope's PSF are needed. Since 
these are not usually available, the scattered light is of- 
ten treated as a non-polarized contribution to the total 
observed light (see Sect. 3). In our study this is partic- 
ularly important because one of the models under study 
{gappy model) postulates the existence of field- free regions 
around the rs = 1 level in the penumbra. These regions 
will naturally produce a non-polarized contribution to the 
total observed Stokes vector. Therefore, there is a po- 
tential risk of not detecting the field-free gaps due to an 
incorrect treatment of the scattered light. 
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Fig. 9. — Azimuthal variation of the magnetic field strength (top), 
line-of-sight velocity (middle) , and inclination of the magnetic field 
(bottom) at rs = 1, for a cut at a radial distance r/Rs = 0.75 
in Fig. 2. The origin of the azimuth angle (abscissa) corresponds 
to the line of symmetry of the sunspot, indicated by the arrow in 
Fig. 1. The scattered light filling factor, Oqs, is also plotted in all 
three panels (dashed lines). 

If our inversions are affected by this degeneracy between 
scattered light and field-free gaps, it is expected that those 
pixels where the intraspines are located show larger values 
for the amount of scattered light retrieved by the inversion 
(ttqs). To study this possibility we have plotted in Figure 9 
the variations of the magnetic field strength (top panel), 
line-of-sight velocity (middle panel) and inclination angle 
(bottom) along an azimuthal cut at r/Rg = 0.75. Other 
azimuthal cuts at different radial distances show very sim- 
ilar behaviors. The values are taken at an optical depth 
of T5 = 1 from the 2-node inversion in Sect. 3.2. This 
plot includes, not only those pixels selected in Fig. 2 as 
intraspines, but all of them. Therefore regions where the 
magnitude of Evershed flow is reduced and the magnetic 
field is more vertical and strong (penumbral spines) are 
also visible. All three panels also show the amount of 
scattered light aqs (dashed lines). There is no particular 
correlation between the location of penumbral intraspines 
(high velocities, weak and very inclined fields) and the re- 
gions where is largest. Similar variations are observed 
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if we plot the values of the magnetic field strength and 
inclination, and line-of-sight velocity, at an optical depth 
of r5 = 10~^. This rules out the possibility that our in- 
versions do not show field-free regions, in the deep photo- 
spheric layers, where intraspines are located at the expense 
of an enhanced scattered light contribution. 

Recently, Orozco Suarez et al. (2007a; 2007b) have pre- 
sented inversions of Stokes spectra measured using Hin- 
ode's spectropolarimeter in the quiet Sun. These authors 
claim that for this instrument it is more appropriate to 
consider a local (unpolarized )scattered light profile. This 
is obtained by averaging the observed Stokes I profiles over 
a small region (about 1 arcsec) around the pixel that is be- 
ing studied. In this different scattered light profile 
is \iscd in the inversion of each pixel. This approach can 
be justified by the fact that the focus of the Narrow-band 
filter (NBI) on Hinode is favored when simultaneous ob- 
servations are carried out with both instruments. In our 
inversions, we have however used a global scattered light 
profile, where we average the Stokes / signal emerging from 
the quiet Sun region far away from the sunspot. In order 
to test whether our results depend on the use of a diff'erent 
scattered light profile, we have repeated our 2-node inver- 
sion using the same approach as Orozco Suarez et al. The 
results are presented in Fig. 10 (cf. Fig. 3). The percent- 
age of intraspinal pixels with dB/dr > is even larger 
than before (90 %). 

Alternatively, Orozco Suarez et al. (2007a; 2007b) point 
out that the most realistic way to account for the scattered 
light would be to consider a local and polarized scattered 
light profile, where not only Stokes / is averaged, but also 
Stokes <5, U and V. We have also tested this possibil- 
ity. Unfortunately, this yields unrealistically high values 
for during the inversion: a^s > 0.9. This indicates 
that the inversion code tries to dominantly reproduce the 
observed Stokes profiles using the scattered light contribu- 
tion. This is not surprising since the neighboring Stokes 
profiles often look very similar to those in the pixel under 
study. Therefore we conclude that this is not a reasonable 
approach when inverting sunspot data. We cannot rule 
out however, that this treatment will work in quiet Sun 
regions. 

As a final test, we have repeated our 2-node inversion 
but neglecting any scattered light: aqs = 0. This test is 
very appropriate because, according to Spruit & Scharmer 
(2006), inversions of spectropolarimetric data fail to de- 
tect field-free regions in the penumbra as a consequence 
of these being already included in the scattered light pro- 
file. If their hypothesis is correct, not accounting for the 
scattered light contribution should uncover these regions 
with _B ~ near rs = 1. Results ignoring the effects of 
the scattered light (Figure 11) are essentially unchanged 
if compared to those where we used a global (Fig. 3) or 
local (Fig. 10) scattered light profile, with 68 % of the 
pixels showing dB/dr > 0. The only difference is the in- 
creased temperatures obtained when we impose aqs = 0. 
In particular we do not see any pixel where the magnetic 
field reaches very small values in the deepest photospheric 
layers. 

Taking into account all tests carried out in this section, 
it seems unlikely that the scattered light can significantly 
bias our magnetic field stratifications, consequently mak- 



ing highly unlikely that we are missing the detection of 
field- free regions near T5 = 1. 

6. CONCLUSIONS 

The uncombed model postulates that penumbral in- 
traspines are characterized by the presence of horizontal 
fiux tubes embedded in a surrounding atmosphere that 
possesses an inclined magnetic field. According to this 
model, looking along these regions should reveal a mag- 
netic field that smoothly decreases at first, but once the 
flux tube contribution starts, the field strength could ei- 
ther increase or decrease. Alternatively, the gappy penum- 
bral model postulates that instraspines correspond instead 
to regions where convectivc field-free gaps penetrate the 
penumbral field. In this case the magnetic field strength 
should also decrease with optical depth at first, but suf- 
fer a much larger drop once the line-of-sight crosses the 
field-free gap. 

In order to differentiate between these two models, we 
have used polarimetric data at very high spatial resolu- 
tion, recorded with the spectropolarimeter on-board of the 
Japanese spacecraft Hinode, to investigate the depth vari- 
ation of the magnetic field strength in the penumbra. We 
have selected a large number (~7500) of pixels that are 
representative of weak and horizontal magnetic fields (i.e.: 
penumbral intraspines) carrying strong Evershed fiows. 
From the inversion of the Stokes profiles at these loca- 
tions we find that the magnetic field strength can either 
increase or decrease with optical depth. A more detailed 
inversion of the average Stokes vector over the selected 
pixels, shows that the magnetic field initially decreases, 
between \0gT5 G [—3, —0.7], but increases thereafter until 
logT5 = 0. 

The gappy penumbral model can explain a smoothly de- 
creasing magnetic field strength only if the T5 = 1 level 
is formed above the field free gap, otherwise a much more 
sudden decrease would be observed as the line-of-sight pen- 
etrates the field-free plasma. A partial solution to this 
problem can be found if we assume that the gap is not 
fully evacuated of magnetic field. However, it offers no 
explanation for about 66 % of the selected pixels, where 
an increasing magnetic field strength with optical depth 
is observed. The absence of field-free gaps, as indicated 
by the inversion, does not in itself imply that there is no 
form of convection present in the penumbra, but rather 
suggests that the convectivc energy transport takes places 
in the presence of a magnetic field (see Zakharov et al. 
2008; Rempel & Schiissler 2008). An example is the roll 
convection proposed by Danielson (1961). 

All inferred stratifications are compatible with the sce- 
nario proposed by the uncombed model. A magnetic field 
that decreases smoothly towards the interior of the Photo- 
sphere can be explained by either a flux tube (of any field 
strength) whose upper boundary layer lies below T5 = 1 
or, if the upper boundary is above T5 = 1, with a mag- 
netic field inside the fiux tube that is weaker than the 
magnetic field a few hundred kilometers above. In addi- 
tion, a magnetic field that increases towards the interior 
of the Photosphere is compatible with a fiux tube with an 
upper boundary layer above rs = 1, but with a stronger 
magnetic field than the one above. This very same con- 
figuration can explain a magnetic field that first decreases 
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and then increases with optical depth, as inferred from the 
inversion of averaged intraspinal profiles. 

We have also studied the effects of the scattered light 
in our inversions. We have seen that any inaccuracies in 
its treatment arc unlikely to be a source of error in the 
stratification of the magnetic field strength. It would be 
very desirable to make a robust confirmation of our find- 
ings, namely dB/dr > in the outer penumbra, for a 
larger number of sunspots at different heliocentric angles 
and including also the disk-ward side of the penumbra. A 
natural extension of this work would be to use the Fe I 
lines at 1.56 /xm (which are formed deeper in the Photo- 
sphere) to confirm the absence of field-free regions around 
T5 = 1. Unfortunately, no such observations exist at the 
spatial resolution needed to resolve pcnumbral iiitraspines 
(~ 0.4"). Indeed, some studies at slightly lower resolution 



(0.6-0.7") have been presented by Cabrera Solana et al. 
(2008), who used simultaneous observations of Fe I 630 
nm and 1.56 fxm recorded with the TIP (Martinez Pillet 
et al. 1999) and POLIS (Schmidt et al. 2003) instru- 
ments. They found that in the outer penumbra, the hor- 
izontal magnetic field component (carrying the Evershed 
flow) was no longer weaker than the more vertical one. 
This can be used as an independent confirmation of our 
work, where we routinely find dB/dr > at large radial 
distances from the center of the sunspot. In addition, flux 
tubes with stronger magnetic field than the one of the en- 
vironment in which they are embedded are also necessary 
to explain certain aspects of the net circular polarization 
observed in the outer penumbra of sunspots (Tritschler et 
al. 2007; Ichimoto et al. 2008). 
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